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Abstract 

Electromagnetic properties of diquarks are investigated in the framework of a color- 
octet Nambu-Jona-Lasinio model, which describes baryons as bound states of di- 
quarks and quarks. We calculate the electromagnetic form factors of scalar and axial 
vector diquark bound states using the gauge-invariant proper-time regularization. 
The nucleon charge radii and magnetic moments are estimated in a simple additive 
diquark-quark model. This picture reproduces the qualitative features such as the 
negative charge radius of the neutron. Within this model axial vector diquarks are 
seen to be important. 
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1. Introduction 

A prime objective of hadronic physics is the understanding of baryons as bound states of 
QCD. Since the direct application of QCD in the nonperturbative regime is extremely dif- 
ficult, models for baryons have been constructed which preserve certain features of QCD 
believed to be essential. Such models include the non-relativistic quark model, in which 
baryons are described as constituent quarks moving in a phenomenological potential, or 
soliton models based on the limit Ac — > oo. In both these approaches the description 
of baryons is essentially of a static nature, which entails serious difficulties due to the 
lack of proper quantization of the translational motion. A fully relativistic description of 
baryons can be accomplished in an approach in which baryons are considered as bound 
states of diquarks and quarks. Evidence for correlated diquark states in baryons is found 
in deep-inelastic lepton scattering |IJ and in hyperon weak decays [Q]. Attempts have 
been made to describe diquarks and baryons in non-local approximations to QCD ||. A 
more tractable approach is based on a local effective quark model, a Nambu-Jona-Lasinio 
(NIL-) model with a color-octet current-current interaction. This model arises as the 
leading term in a current expansion of the quark-quark interaction and can also be mo- 
tivated semiclassically Q. It describes the spontaneous breaking of chiral symmetry, i.e., 
the dynamical generation of quark masses. Furthermore, after Fierz transformation of the 
interaction the color-singlet part of the model can be bosonized exactly and provides a 
realistic description of meson properties H, ||, 0]. Diquark bound states have been consid- 
ered in |8], H|. In order to study baryons, the NJL model with color-octet interaction has 
been transformed into an effective theory of meson and baryon fields by use of functional 
integral techniques |I0|. In this theory bound states of diquarks and quarks occur due to 



quark exchange. In |Tl[ the Faddeev-type equation for the binding of a scalar diquark 
and a quark has been solved using a static approximation for the quark exchange inter- 
action. The resulting masses of the | + -baryon octet are in satisfactory agreement with 
experimentS. 

Our aim is to extend this approach to the study of the low-energy electromagnetic 
structure of baryons. As a first step towards describing baryon form factors in the effective 
hadron theory derived from the NJL model, we calculate here the on-shell electromagnetic 
form factors of the constituent diquarks. The relevant diquark bound states occur in the 
scalar (0 + ) and axial vector (1 + ) channel. The diquark form factors are a crucial ingredient 
in the baryon form factor and contain information about the sizes of the correlated diquark 
states. To ensure gauge invariance, the NJL quark loop is regularized using Schwinger's 



proper-time method fl3fl . We also compare with the results obtained with the usual sharp- 
cutoff regularization. With the diquark form factors we then estimate the nucleon charge 
radii and magnetic moments in an additive model of extended diquarks and pointlike 
valence quarks, using the spin-flavor wave functions of the non-relativistic quark model. 
Such a picture may be regarded as a crude approximation to the baryon wave function. 

2. The model 

The basis of our description is an effective quark Lagrangian of the NJL type, 

C q = m -m°)q-ig^f A , (1) 
where = q'j^^X A q is the color-octet current of the quark field, g an effective coupling 



a Recently, direct numerical solutions to the Faddeev equation for the nucleon have been obtained [|Tj 
which support the static approximation pj|. 
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constant and m° the current quark mass matrix. By way of Fierz transformations in 
color, flavor and Dirac space the interaction in eq.(|l|) can be rewritten equivalently as 
an interaction in the color-singlet (lc) quark-ant iquark and color-antitriplet (3c) quark- 



quark channel [10], 



-\gjtf A = 9i(qA a q)(qA a q) + g 2 (qT a Cq T )(q T CT a q). (2) 

Here, C = i'j 2 '-/ is the charge conjugation matrix and the vertices A a , T a are matrices in 
color, flavor and Dirac space, 

a- = 10 (£)<>•, r=(*)$)o-, -14V), (3) 



where 



and (e A )sc — £abc are the antisymmetric generators of the 3-representation of SU(3)c- 
From the Fierz transformation one obtains g\ = gi = g. However, the two terms in eq.(0) 
are independently chirally invariant, so that chiral symmetry does not preclude different 
effective coupling constants in the meson and diquark channel. 

By introducing collective meson and baryon fields, the generating functional defined 
by the quark lagrangian, eq.([I|), has been rewritten as an effective hadron theory |TIJ. In 



the course of this reformulation one also introduces collective fields for diquarks, which 
form building blocks for baryons. One obtains an effective action of meson and diquark 
fields, (p a and A a , coupling to a baryon source, x> 



S = ~J d 4 x (0> a + A* Q A Q ) - \i Tr A log G- 1 - J d 4 x ( X X T ) Ob [A] ^ j 



.(5) 



Here, G is the quark Green's function in the Nambu-Gorkov formalism of super conduct i- 
vit\ 



G 1 = ( _(^it j > G 1 = i $ -m -(/)-- I (fx XX, (6) 



and the collective meson and diquark fields couple to the quarks through the vertices 
= (p a A a and A = A Q r a , A = A*°T a . In particular, the vacuum value of the scalar 
meson field generates the constituent quark mass, M. The form of the baryon propagator, 
Gb[A] is given in [|1^]. After integrating out the diquark fields one obtains from eq.(|5]) an 
effective hadron theory, in which diquarks and quarks interact through quark exchange 
IIJ and form bound baryon states . 



The quark loop in eq.(|5]) is defined with a cutoff procedure, which is an important 
physical ingredient of this effective model. A simple regularization scheme is a sharp 
Euclidean cutoff, applied to the momentum integrals resulting from the expansion of 
Tr logG -1 in meson and diquark fields. This method in general breaks gauge invariance, 



b Following ref. |hJ we denote 2 x 2-matrices in the Nambu-Gorkov formalism by bold sans-serif 
letters. 
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when an electromagnetic field is introduced. A gauge invariant regularization is provided 



by the proper-time method of Schwinger [13 , which defines the real part of the quark 



loop, after continuation of G to Euclidean space, as 

roc da , 

Re Tr A log G E = \ \ — Tr exp(-s G^G S ). (7) 

J A' 2 S 

The imaginary part of the Euclidean quark loop, which generates anomalous (intrinsic 
parity-violating) meson and diquark processes, is finite and left unregularized ||. 

3. Diquark electromagnetic form factors 

An important ingredient in the electromagnetic structure of baryons are the electromag- 
netic form factors of the constituent diquarks. To describe the electromagnetic interactions 
of diquarks, we couple an electromagnetic field to the quark fields of eq.(|J), by way of 
minimal substitution, — > i(j) — QA^, where Q = |e(A 3 + ^A 8 ) is the quark charge 
matrix. We then expand the effective action, eq.([5|), in the background of the electro- 
magnetic field to second order in the diquark field. (Here, the baryon source is set to 
zero, x = 0, an d the meson field is left at its vacuum value.) In Minkowski space and in 
momentum representation, this expansion is of the form 



(2k) 

+ /(04 j^Ki-V-\q)^ip-\q)AM^^<ll (8) 

where T> is the diquark propagator and T the electromagnetic vertex function. The masses 
of the diquark bound states are determined as the zeros of the inverse propagator. On-shell 
diquark form factors are then obtained by evaluating the vertex function for appropriate 
incoming and outgoing four-momenta p±\q on the diquark mass shell and normalizing the 
fields in eq. (H) to unit residue of the propagator. Note that the electromagnetic couplings 
of the diquark fields come entirely from the quark loop of eq. (§) . The diquark propagator 
and electromagnetic vertex function are diagonal in color and of the formi 



(9) 



v a \ / \b s 



^( P , q ) = ^E(v) (t) (Qi4 b *M-Qi#*(-p>9))- ( 10 ) 



v 2 / .. V 2 , 

In cutoff regularization, the functions if? and jf^ are given by the loop integrals 

ifiP) = ¥f A2 ^^^[G i (k-lp)0 a G j (k+lp)0\ (11) 

J?^(P, q) = \i / A2 (04 tr Di- [Gi(k - \p)0*G 3 (k + \p - \q)^G 3 (k + \p+ \q)0\ 

where Gi(p) = (j> — Mj) -1 is the quark propagator. In proper-time regularization cor- 
responding expressions are derived using standard techniques for the expansion of the 
time-ordered exponential, eq.(^) Only the final expressions will be given below. 



c The symbol g 3 * means the metric tensor if a, b are Lorentz indices, otherwise 5^ . 
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Table 1: The masses of the scalar ud- and the axial vector uu-diquark, m + and 
mi+, and their electromagnetic charge radii, for various values of the effective diquark 
coupling constant, g 2 . Also shown is the magnetic moment of the axial uu-diquark, 
Hi+ , in units of 2Q u e/2M. Proper-time regularization is used; the parameters are 
M = 400 MeV, A = 630 MeV, m° = 17 MeV. 



We now evaluate eqs.(|9], |H]) for the scalar and the axial vector diquark channel, which 
are relevant to the description of baryons. We use throughout the proper-time regulariza- 
tion and only compare the numerical results with those obtained with cutoff regularization. 
Since we shall consider only the nucleon later we assume the isospin limit, M u = Ma = M. 
Let us first consider the scalar Mo?-diquark (O a = 275). In proper-time regularization one 
finds, after continuation to Minkowski spaceB, 

I 0+ (p) = p 2 A (p 2 )+2M 2 B, (12) 

with 

- i^7, "H"' '" "t*~ a * )- (13) 

1 / M 2 \ 

B = 

The second term in eq.(|I~2"D is related to the quark condensate, (uu) = (dd) = 4NqM 3 B. 
The scalar diquark mass, m +, is shown in table 1 for different values of the effective 
diquark coupling constant, g 2 - Here, the coupling in the meson channel, gi, is determined 
by the constituent quark mass through the gap equation, the cutoff is fixed by fitting the 
pion decay constant, and m° is determined from the pion mass. Diquark couplings of 
92/ gi ~ 2 are required to reproduce the masses of the spin-| baryons if quark exchange 
is included JTTJ. A ratio ^/^i ~ 2.5 is needed in an additive diquark-quark model, see 




below. Similar values are also required to obtain sufficiently bound axial diquarks. 

To obtain the electromagnetic form factor for an on-shell scalar diquark of mass m + , 
we evaluate the vertex function, eq.(10), for incoming and outgoing momenta p± with 
(p ± |g) 2 = mg + , which entails p ■ q = 0. On the mass shell, the vertex function is 
transverse to the photon momentum, 

J^ + (p,q) = 2F 0+ ( q 2 )p», (14) 

and the normalized diquark form factor is defined as fo+(q 2 ) = F +(q 2 )/Z +, with Z t 
3/ dp 2 J + |„2 =m 2 . In proper-time regularization one finds 



0+ 



"0+ 



F 0+ (q 2 ) = A (m 2 0+ ) + (m 2 0+ - h 2 )^ 2 ) + WC 2 (q 2 ) (15) 



d We suppress the flavor indices on I 3 ^ and J 3 ' 3 in the following. 
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Figure 1: The electromagnetic charge form factor of the scalar ud-diquark, fo+{q 2 ), 
(solid line) and the axial uu-diquark fl+(q 2 ) (dashed line), for m + = 318 MeV and 
mi+ = 764 MeV fa/gi = 2.5, M = 400 MeV). Shown are the results for proper-time 
regularization. The curve for the pion form factor would nearly coincide with the one for 
the scalar diquark. 



where 

c,M 2 ) = i^i Sl dp Io~" daXl ' 2 exp( ~y 2 VA2) , (16) 

X 1 = a, X 2 = 1, Y 2 = M 2 — a(l — a)ml+ — (3(1 — a — P)q 2 

The scalar diquark form factor as a function of the photon momentum, q 2 , is shown in fig. 
1, for a diquark mass of m + = 318 MeV. In particular, /o+(0) = 1, i.e., the total charge 
is conserved as a consequence of the gauge-invariant regularization. The scalar diquark 
r.m.s. charge radius, (r 2 ) + = — 6 d/dq 2 f + \ q 2=o is given in table 1. 

In sharp cutoff regularization the scalar diquark mass and form factor are evaluated 
by continuing the integrals (O ) to Euclidean space, introducing Feynman parameters 
and applying the cutoff after a shift of the integration variable. Due to the breaking of 
gauge invariance by the cutoff, /o+(0) ^ 1 if "m\+ > 0, so that we have divided fo+(q 2 ) by 
/o+(0) to obtain meaningful results. Determining the parameters as above in the meson 
sector, somewhat larger ratios gijg\ are required to obtain the same diquark masses as in 
proper-time regularization. Nevertheless, comparing the charge radii for scalar diquarks 
of the same mass, the results are quite close in both regularization schemes, cf. table 2. 
The virtue of the gauge-invariant proper-time regularization, eq.(|7|), is that it is applied 
at the level of the action, rather than at the level of individual loop integrals, which is 
necessary for the electromagnetic Ward identities to be satisfied. In particular, with this 
method the electromagnetic vertices are also defined off the diquark mass shell. 

It is instructive to compare the form factor of the scalar wd-diquark to that of the 
charged pion [ 14]| . The only difference between the diquark and the meson vertex function 



is the sign of the quark charges in eq.flTOD and a factor Nq multiplying the quark loop. 
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proper-time 


sharp cutoff 


m + /MeV 


mi+ /MeV 


(r 2 )^ 2 /fm 


(r 2 )\{ 2 /fm 


mi+ /MeV 


(r 2 )^ 2 /fm 


(r^f/fm 


400 


778 


.48 


.90 


785 


.43 


1.02 


300 


761 


.48 


.74 


770 


.42 


0.79 



Table 2: Comparison of the diquark charge radii for proper-time and cutoff regularization. 
Here, the ratio g^l g\ is chosen independently in both schemes in order to obtain the given 
scalar diquark masses. 



However, since the same factor iVc occurs also in the meson propagator, this factor cancels 
if we consider the normalized on-shell pion form factor. Thus, in the isospin limit the 
pion form factor is given by the expression for the diquark form factor, eq.(15), evaluated 



with m,Q + — > m 2 . The curve for the pion form factor as a function of q 2 is almost identical 
to that of the scalar diquark shown in fig. 1. 

One may consider the scalar diquark (or pion) r.m.s. charge radius in table 1 as a 
function of the bound state mass. For strongly bound diquarks the charge radius is 
essentially independent of the diquark mass and thus practically identical to that of the 
pion. It is therefore, to good approximation, given by the result of a gradient expansion 
around p 2 = |T5| , which in proper-time regularization reads 

(r 2 )o + = ^exp(-M 2 /A 2 ). (17) 

In the other extreme, near m = 2M, where the diquark would become unbound, the 
charge radius grows like (r 2 ) + ~ (m — 2M)~ X . This is the behaviour expected for a 



weakly bound state below a continuum threshold ||16|| . Such divergent behaviour would 
be absent in a model which incorporates quark confinement. 

Let us now consider the axial vector diquark (O a ' b = ry p '°" / 'y/2) . With the gauge- 
invariant proper-time regularization and in the absence of flavor symmetry breaking, the 
axial diquark polarization operator deriving from the quark loop is transverse, 

I£(p) = p 2 A 1 (p 2 )(g^-^), (18) 

with 

„ / 2n 1 f 1 , / s^f M 2 -a(l-a)p 2 \ , , 

Mp 2 ) = — j o da2a(l-a)V[o, A_ IE. j . (19) 

The longitudinal part of eq.(§) is thus given entirely by (3/2g2)g pa . Axial diquark masses 
are given in table 1 for various values of g^. In contrast to the scalar diquark the axial 
diquark is only weakly bound for reasonable values of g^. The electromagnetic vertex func- 
tion, eq. (|TUD , for on-shell axial vector diquarks of incoming and outgoing four-momenta 
p ± \q, with (p ± \q) 2 = m\ +) is of the form JHJ 

J&faq) = 2gPy^F l \{q 2 )-[g^{p+\qr + g^{p-\qY}F^{q 2 ). (20) 

In proper-time regularization one findsi 

Ft + {q 2 ) = ^(q^-H^q 2 ), (21) 



e In eq.(|20|) we have omitted a term oc (p— ^q) p (p+ \q) a p p , since we are interested only in the diquark 
charge and magnetic moment form factors. 
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F? + (q 2 ) = |AoK + )-iA (g 2 )-ig 2 C 1 (g 2 ) + X + C '2(g 2 )-^ 2 (g 2 ), (22) 

where 

#1,2 = ^ £daJ^ a df3W h2 r , W 1 = a, W 2 = 2(3, (23) 

and F 2 and F +(g 2 ) are as defined eqs.([U], |15) and evaluated with m + — > mi+. The 
normalized electric and magnetic form factors of the axial vector diquark are obtained as 
f^(q 2 ) = F e ^(q 2 )/Z 1+ , where Z x + = d/dp 2 (p 2 A 1 )\ p 2 =m 2 + . In particular, ff+(Q) = 1 and 
f^+(0) is the magnetic moment of the zj-axial vector diquark in units of (Qi + Qj)e/2mi+ 
(i,j = u,d). The electric form factor describes the charge distribution in the axial di- 
quark and is shown in fig. 1. The charge radius of the axial diquark is larger than that 
of the scalar diquark, see table 1. This fact is mainly due to the weaker binding of the 
axial diquark. However, even strongly bound axial diquarks would be larger than scalar 
diquarks of the same mass. In gradient expansion at m + = = one would find 
(r 2 )x+ = 7^(r 2 )o+. Larger charge radii for the axial diquarks are also found in cutoff 
regularization, cf. table 2. The axial diquark magnetic moment is shown in table 1, in 
units of the sum of the constituent quark magnetic moments, (Qi + Qj)e/2M. The di- 
quark magnetic moment is seen to be lowered compared to the non-relativistic value as a 
consequence of the diquark binding. 

4- Estimating the baryon form factor 

Given the diquark electromagnetic form factors we now wish to describe the electromag- 
netic interactions of baryons as bound states of constituent diquarks and quarks. To 
calculate the baryon electromagnetic form factor in the effective hadron theory [l(J based 



on the NJL model, one would extract from eq.(|5|) the baryon electromagnetic vertex, 
which describes the coupling of the photon to the diquarks, the quarks and to the quark 
exchange. One would then take the matrix element of this vertex between on-shell baryon 



states, i.e., solutions of the baryon Faddeev equation JTT]. Since a full calculation of the 
form factor using bound state wave functions is rather involved, it is worthwhile to first 
explore a simpler picture. Our intention here is to see to what extent the nucleon charge 
radii and magnetic moments can be described in an additive diquark-quark model ||. 
This approach has been successful in the description of baryon masses. It is based on 
the spin-flavor wave functions of the non-relativistic quark model, in which two quarks 
are coupled to a scalar or axial diquark. The diquarks are then regarded as on-shell 
constituents, whose intrinsic properties such CIS db IXlclSS, db charge radius or a magnetic 
moment, are derived from the NJL model. Note that as a consequence of the Pauli prin- 
ciple the scalar diquark is antisymmetric in flavor (ud) and has isospin 0, while the axial 
diquarks are symmetric (uu,ud,dd) and have isospin 1. 

We first consider the nucleon electromagnetic charge radii. In a simple approximation 
one may neglect the orbital motion of diquarks and quarks, i.e., take the orbital wave 
functions of diquarks and quarks in the baryon Breit frame to be 5-functions centered at 
the baryon center-of-mass. The charge density inside the baryon is then the sum of the 
extended diquark charge distributions and the pointlike valence quark charge densities, 
which do not contribute to the charge radius. From the quark model wave functions one 
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then obtains the relatione 

(r\ = |(|(r 2 )o + + (r 2 )i + ), (r 2 )„ = i(|(r 2 ) 0+ -|(r 2 )i + ). (24) 

Note that the diquark charge radii have been defined as derivatives of the electric form 
factor normalized to 1 at q 2 = 0, so that the diquark charges are absorbed in the numerical 
factors in eq.fl2"3|). We now choose the diquark coupling, g 2 , such as to fit the proton mass, 
using the formula of the diquark-quark model, m P)n = |(m + +mi+) +M |J. From table 1 
we see that this requires a value of g2/gi = 2.5. With the diquark charge radii from table 1 
eq.fl2"4|) gives (r 

2) = (0.57 fm) 2 and (r 2 ) n = -0.06 fm 2 , which is in qualitative agreement 



with the experimental values of (r 2 ) pcxp = (0.85 fm) 2 and (r 2 ) njexp = — O.llfm 2 . This 
simple picture naturally explains the negative charge radius of the neutron by the fact 
that the axial diquark charge radius is larger than the scalar one. Note that this qualitative 
result is independent of the diquark masses. 

If one included only scalar diquarks in the above estimate, one would obtain instead 
of eq.(^) (r 2 ) p = (r 2 ) n = |(r 2 ) +, i.e., the neutron charge radius would come out equal 
to the proton one. This underscores the importance of the axial vector diquarks in the 
additive diquark-quark model. We remark, however, that in a model with only scalar 
diquarks a negative charge radius for the neutron can be obtained if the orbital motion 
of diquarks and valence quarks is taken into account, as has been demonstrated in the 



framework of the non-relativistic quark model [18|. The question of the relative impor- 
tance of the intrinsic size of the axial diquarks versus the effects of the diquark-quark 
orbital motion can only be answered a posteriori from the full relativistic bound state 
wave function containing both scalar and axial diquark components. Nevertheless, it is 
encouraging that a simple additive picture with axial diquark constituents of a size close 
to that of the proton is capable of reproducing both the nucleon mass and the charge 
radii reasonably well. Moreover, in the non-relativistic quark model it can be seen that 
the estimate of the proton charge radius in eq. ([24]) is increased by the diquark-quark 
relative motion, while the neutron charge radius remains essentially unchanged unless 
S'?7(6)-breaking in the spin-flavor wave function is taken into account [|T^]. This strongly 
suggests that the estimates of the additive model will be improved by a full calculation 
using Faddeev wave functions. 

Finally, we want to discuss the nucleon magnetic moments in the diquark-quark pic- 
ture. Since magnetic moments are described very well by the non-relativistic quark model, 
it is important to see how far the predictions of the diquark-quark picture differ from those 
of the quark model. Rewriting the quark model spin-flavor wave functions in terms of 
diquarks and quarks one obtains 

= 7^ (| + + §/^o+-i+) , Mn = ^ (~l ~ 9^1+ _ |Mo+-i+) • (25) 

The first term here is the valence quark contribution. Furthermore, is the magnetic 
moment of the (ij)-axial diquark (i, j = u, d) in units of (Qi + Qj)e/2M, and ^to+-i+ is the 
transition moment from the u<i-scalar to the lid-axial diquark in units of (Q u — Qd)e/2M, 

((ud) +,S z = 0\fi lz + fj, 2 z\{ud) 1+ ,S z = 0) = /i +-i+(<5« - Q d )(e/2M). (26) 



Here, the isospin limit is assumed, so that the axial uu, ud- and dd-diquarks all have the same charge 

>!+• 
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Here, /ii z , [L 2z are the z-components of the quark magnetic moment operators. In the case 
of SU(6) symmetry one has = yU +-i+ = 1, which leads to the well-known results 
/i p = e/2M, /i n = — |e/2M. We now consider eq. fl25|) in the sense of the additive diquark- 
quark model outlined above. The magnetic moment of the axial diquark, fi^, has been 
calculated above and is given in table 1. Transitions between scalar and axial vector 
diquarks due to the electromagnetic field are also found with the effective diquark action 
derived from the NJL model, eq.([5|). This is an anomalous (intrinsic parity-violating) 
process, which is generated by the imaginary part of the quark loop; its meson analogue 



is the process p — > ir^y. Since the definition of /z +-i+ in eq. fl26D involves vanishing photon 
momentum, it is not possible to describe /io+-i+ as a transition between on-shell scalar 
and axial diquarks of unequal mass. To obtain a rough estimate one may evaluate this 
vertex for a scalar diquark mass equal to that of the axial diquark. In this case the 
transition part of the vertex function, eq.(|i~0"l), is given by 

J^-i+Cp,?) = ^e^VM-^o^i^p 2 ), (27) 
i , 2M 2 



Fa+~i+{p 2 ) = f da 
Jo 



o M 2 -a(l-a)p 2 ' 

and the transition moment corresponds to // +-i+ = F +_ 1 +(m 2 + )/V2Z +.Z 1 +. For the 
axial diquark masses considered in table 1 /xq+-i+ is close to its non-relativistic value of 
1; yUo+-i+ = 0.97 for g 2 / ' g\ = 2.5. For this value of g 2 /gi = 2.5 we obtain from eq. (|25"D 
fi p = 0.96 e/2M, \i n = — 0.65e/2M. Note that the constituent quark masses used in 
the NJL model, M ~ 400 MeV, are larger than those of the quark model, M = |m p ~ 
300 Me V, so that come out roughly 3/4 smaller than the quark model prediction. 
We emphasize, however, that a proper treatment of the transition contribution would 
require the full baryon wave function with off-shell diquarks. Note that if one included 
only scalar diquarks in the additive diquark-quark model, the magnetic moments would 
be \i p = ~e/2M, p n = — |e/2M. The same result would be found if one included 
axial diquarks but neglected axial-scalar transitions. These estimates again show the 
importance of axial diquarks in the additive model. 

5. Conclusion 

We have calculated the electromagnetic form factors of scalar and axial vector diquark 
bound states in a color-octet NJL model. For realistic diquark masses the scalar diquark 
charge radius is close to that of the pion. The axial diquark is weakly bound; its size is of 
the order of that of the proton. The nucleon charge radii can be described qualitatively 
in an additive diquark-quark model, which may be viewed as a crude approximation to 
the baryon wave function. The inclusion of axial vector diquarks in this scheme is seen 
to be crucial. For the nucleon magnetic moments scalar-axial diquark transitions have to 
be taken into account at least approximately. It remains to be seen to what extent the 
results of the diquark-quark model are improved by a dynamical calculation using the 
relativistic baryon wave functions from the Faddeev equation. 
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